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We discuss the influence of two-dimensional hexatic order on capillary waves and undu- 
lation modes in spherical and cylindrical geometries. In planar geometries, extended bond- 
orientational order has only a minor effect on the fluctuations of liquid surfaces or lipid 
bilayers. However, in curved geometries, the long-wavelength spectrum of these ripples is 
altered. We calculate this frequency shift and discuss applications to spherical vesicles, liq- 
uid metal droplets, bubbles and cylindrical jets coated with surface-active molecules, and 
to multielectron bubbles in liquid helium at low temperatures. Hexatic order also leads to a 
shift in the threshold for the fission instability of charged droplets and bubbles, and for the 
Plateau-Rayleigh instability of liquid jets. 
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I. INTRODUCTION 

In two dimensions the melting from a crystal to an isotropic liquid can be a two-stage process 
|Q], driven by the sequential unbinding of dislocations JHQl and disclinations [1]. At low temper- 
atures T < T m dislocations are suppressed due to their cost in elastic energy. However, their free 
energy decreases with increasing temperature. At a temperature T = T m , the quasi-long-ranged 
translational order of the crystal is destroyed by the dissociation of dislocation pairs. This transi- 
tion leads to an intervening hexatic phase, which still exhibits extended orientational correlations. 
The unbinding of disclination pairs sets in at a higher temperature T = Ti. In this second transition 
the quasi-long-ranged orientational order of the hexatic phase is destroyed, leading to an isotropic 
liquid. This mechanism allows the melting transition to be continuous in contrast to the first-order 

n 

melting (directly to an isotropic fluid) predicted by Landau [3]. 

Several experimental systems have illuminated the nature of two-dimensional melting. A nearly 
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ideal system is electrons on helium ||4j, |5l |6J . The electrons are trapped on the surface of liquid 
helium by a submerged, positively charged capacitor plate. Their separations are rarely less than 
1000A, so the in-plane physics is that of classical particles with a repulsive 1/r potential. The 
liquid helium does not freeze at low temperatures, so it is possible to cool the electrons on this 
liquid substrate well below a sharply defined 2d freezing temperature (T m ~ 0.5K) On the 
theoretical side, important parameters such as the 2d shear modulus and dislocation core energy 

R 

are easily calculated with this potential [7]. Computer simulations [8] reveal a shear modulus 
which appears to drop to zero at the melting temperature. Measurements of the shear modulus 
|E] and specific heat [5] are consistent with a continuous dislocation mediated melting transition. 
However, in these experiments it is difficult to determine experimentally if hexatic order and a 
disclination unbinding transition are in fact present above T m . 

However, experimental evidence for hexatic order has been found in a variety of other systems, 



including free standing liquid crystal films and Langmuir-Blodgett surfactant monolayers IjlCjl. 
A hexatic-to-liquid transition has been observed in two-dimensional magnetic bubble arrays 111 ill . 
Furthermore, there is strong support for two-stage continuous melting from recent experiments on 
two-dimensional colloidal crystals ifl^ fisl. In this case, the colloids can be directly imaged by 
video-microscopy, thus allowing a precise test of the theory. 

The modest time scales available even on the fastest computers make equilibration in Monte 
Carlo or molecular dynamics simulations of two-dimensional melting difficult. However, there 
is now evidence via computer simulations for continuous melting and a narrow sliver of hexatic 
phase for hard disks [14] and for particles interacting with a repulsive 1/r 12 potential Jl^ll . There 
are also indications of defect mediated melting transitions for the familiar Lennard- Jones 6-12 pair 

potential y. 

In the above experiments, order was typically probed via diffraction or by direct measurement 
of correlation functions in real space. It is difficult to use these methods, however, when hexatic 
order is present on a curved surface, such as a sphere or a cylinder. Examples where hexatic and 
crystalline order might be present on a sphere include "liposomes", i.e., closed vesicles composed 
of lipid bilayers JlJ, the surface of liquid metal droplets confined in Paul traps Jl^j . and 
multielectron bubbles submerged in liquid helium Il20ll . 

Hexatic order in spherical liposomes seems likely because flat two-dimensional planar layers of 
lipids such as DMPC (Dimyristoyl phosphatidylcholine) [2lj], similar to free standing liquid crys- 



tal films [22], can exist in a variety of states with different degrees of positional and orientational 
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order. Examples include fluid, smectic C, hexatic and crystalline phases. 



Curved hexatic order may also arise on droplets. Celestini et al. 112311 have found evidence 



from computer simulations for extended orientational correlations at the flat surface of supercooled 
heavy noble liquid metals, such as Au, Pt, or Ir. These metals have a general tendency to reduce the 
interatomic distance at the surface. Upon supercooling this effect is enhanced, leading, generally 
to a two-dimensional crystalline surface layer. Under suitable cooling conditions, hexatic order 
can also appear on the free surface of undercooled liquid metals [23]. Surface hexatic order may 
also occur on water droplets coated with surface-active molecules, given that there are already 
extensive observations of this type of order in Langmuir-Blodgett monolayers lioll . 

Finally, one might expect spherical hexatic and crystalline order in multielectron bubbles in 
helium. These arise when the helium surface undergoes an electrohydrodynamic instability at 
high electron densities. The surface then develops a regular array of dimples, each containing 10 6 
electrons or more. As the electric field increases these dimples deepen until eventually electrons 
break through the interface. After subduction, large numbers of electrons (10 5 — 10 7 ) then coat the 
inside wall of a large ( 1 0-1 00/im radius) sphere of helium vapor. These multielectron bubbles have 
been observed to move through the helium after their creation above a protruding anode J20I 
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2511 . They are stable at low electron densities since the Coulomb pressure can be compensated by 



the surface tension of the helium liquid- vapor interface. However, if the electron density becomes 
too high the electrostatic repulsion exceeds the balancing force and the bubble undergoes fission 

HQ. 

One might hope that bond-orientational order in a flat membrane or interface could be detected 
by its effects on the dynamics of undulation modes or capillary waves. Unfortunately, hexatic order 
couples only to the Gaussian curvature W28I1 . which vanishes for a simple sine wave deformation of 
a flat membrane or an interface (cf. Eqs. (l27t and (l84l below). The situation is different, however, 
when these excitations are superimposed on a nontrivial background geometry such as that of a 
sphere or a cylinder. In a recent short communication [29] we have determined the effect of hexatic 
order on the undulation modes and capillary wave excitations for the spherical systems described 
above. The frequency shift is large for liposomes with hexatic order. Observable effects could 
also occur for liquid metal droplets, surfactant coated water drops and in multielectron bubbles in 
helium. In this paper we describe hexatic dynamics on spherical surfaces in detail and extend the 
theory to include cylindrical geometries. 

Cylindrical geometries could be realized by, e.g., coating a liquid jet with a hexatic monolayer. 
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These jets will (similar to conventional liquids) undergo the well-known Plateau-Rayleigh shape 
instability as soon as their length reaches a critical size. However, the stiffness associated with 
extended orientational correlations shifts the threshold of this instability and alters the decay of 
the cylinder into a chain of droplets. Cylinders provide also an example where hexatic order is 
perfectly compatible with the underlying geometry: since the Gaussian curvature of the cylinder 
vanishes, no disclination defects are present in the ground state to complicate the analysis. 

The remainder of this paper is organized as follows. First, we analyze the influence of surface 
hexatic order on spheres by considering liquid droplets (Sect. HTJ). Then, in Sect, [in] we a Pply our 
analysis to cylindrical geometries. Here, we concentrate on liquid jets and we determine the effect 
of hexatic order on the Plateau-Rayleigh instability. Next, we compute the shifted instability due 
to hexatic order in multielectron bubbles (Sect. llVl) . As a last application we investigate in Sect.lVl 
the influence of hexatic order on the undulations of spherical vesicles. Finally, in Sect. |VH we 
elaborate in detail on experimental consequences of our work. 



II. DYNAMICS OF LIQUID DROPLETS WITH HEXATIC ORDER 

We first discuss liquid droplets with surface hexatic order. The equilibrium shape minimizes a 
droplet free energy Fd given by contributions from an interfacial energy and the hexatic degrees 
of freedom 

F^Fi + Fh^a J dA+ l -K A J dA D^D'nj, (1) 

where we use the summation convention throughout and a denotes the surface tension of the 
interface of the liquid droplet. For a general manifold with internal coordinates x = (x l ,x 2 ), 
the surface element is given by dA = a/ g(x)d 2 x, where g(x) is the determinant of the metric 
tensor gij(x). For an undeformed sphere with radius R , x = (9, <p) with polar coordinates 6 and 
if and dA = R$ sin 9 d6 dip. The quantity n is a unit vector in the tangent plane with n^n* = 1 



which identifies (modulo 27r/6) the long-range correlations in the hexatic bond directions [28]. 
Here, D^nP = g^ k Dink, where g^ is the inverse of gij. The operator Di denotes a covariant 
derivative with respect to the metric o^-. Thus, = div? + T J ki n k , where the T 3 ki are Christoffel 



symbols of second kind. See, e.g., IB0I1 . Close to the melting temperature T m , the hexatic stiffness 



K A ~ / clq) 2 , where £t is the translational correlation length, ao is the particle spacing and 
E c is the dislocation core energy yj]. The ratio KaI^bT jumps from an universal value 72/n to 
zero when the hexatic melts into an isotropic liquid at T = T-~ (see Fig.[TJ Jjl. 



5 




FIG. 1 : The hexatic stiffness Ka as function of the temperature T. Ka diverges near T m and jumps 



discontinuously to zero at T with lim T _™- Ka(T) = llhsTji: 



31]. 



Droplets have a nearly constant volume V and the corresponding constraint could be included 
in the free energy However, here it is easier to first consider shape fluctuations which explicitly 
keep the volume fixed. Finally, in Eq. dU) we neglect effects arising from gravity since we only 
consider droplets with radii R <C l c much smaller than the capillary length l c , which is of order 
millimeters or more for typical droplets in the earth's gravitational field. 

In the following we will consider the dynamics of shape fluctuations of spherical droplets. As 
the droplet deforms its geometrical properties change, i.e., the metric and the mean and Gaussian 
curvature are altered. The free energy dU) thus has a functional dependence on the underlying 
droplet shape. This dependence can be treated most efficiently by parameterizing the surface of 
the droplet by its surface vector R. For a sphere, we have 

Rfa 1 , x 2 ) = R(9, if) = Rq(9, <£>) (sin 9 cosy?, sin^sinyj, cos#) . (2) 

See Appendix lAl for important geometrical quantities such as the metric tensor gij(x), the mean 
curvature H(x) and the Gaussian curvature G(x) in terms of the surface vector R^x 1 , x 2 ). 

The fundamental assumption which underlies the hexatic free energy discussed above is that 
the configuration of minimal elastic energy corresponds to a vector field n where n (x + dx) can 
be obtained from uq(x) by parallel transport of n . On a sphere however, curvature introduces 
"frustration" since parallel transport of n along closed loops on the surface leads to a rotation of 
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n. Because of this frustration the ground state of hexatic order on a sphere has at least 12 positive 



disclinations. This constraint is a consequence of the Poincare index theorem 13211 . which states 
that a vector field on a surface with genus g and Euler characteristic E = 2(1 — g) must have 
singularities with total vorticity 2ttE. As a consequence, order which is identified by a vector 
order-parameter field on a curved geometry frustrated by a nonzero integrated Gaussian curvature 

G = J dAG(x), (3) 

always has topological defects Il33t1 . On a sphere with g = and E = 2 hexatics must have a 
minimum of 12 defects with charges 2tt/6 W34I1 . 

The energy of an isolated disclination in a hexatic diverges logarithmically in flat space. How- 
ever, this energy is reduced due to screening by the Gaussian curvature of the sphere. This point 
can be made more precise by introducing a local bond-angle field 6, the angle between ft and some 
local reference frame. The transverse part of 9 is then connected with the disclination density 
HQ. As shown in Appendix HO the elastic free energy associated with hexatic order can then 



be written as 135 



F h = ~K A j dAj dA \G{x) - s(x)} T(x, x') [G(x') - s(x')} . (4) 

Here, T(x, x') is the inverse Laplacian on the sphere (see Appendix El), and s(x) the disclination 
density fl33fl - 

1 N d 

^==^2qi$(x - (5) 



six 



x 



1=1 



with Nd disclinations of charge g$ = ±2tv/6 at positions X{. In deriving Eq. © we have assumed 
that the regular part of the bond-angle field 9 reg relaxes rapidly on the time scale of shape deforma- 
tions. In Sect. |V] we will show that this assumption is indeed justified for the systems considered 
here. Finally, it should be emphasized that Eqs. © and © hold for arbitrary geometries, not just 
for that of the sphere (cf. Appendix |Bl). 

The defects minimize F h by arranging themselves to approximately match the Gaussian curva- 
ture, which is G(x) = 1/Rq for a rigid sphere of radius Ro- Deep in a hexatic phase on a sphere, 
we expect Nd = 12, corresponding to 12 fivefold disclinations which lie on the vertices of an 
icosahedron. With polar coordinates such that there are disclinations at the north and south pole, 
the 12 defect locations entering Eq. © are given by 

f 2irk 7r 2j[k 1 

(9k, fk) e < (0, 0), (7, ^r- )o<fc<4, (?r - 7, - + )o<fc<4, (tt, 0) > , (6) 
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FIG. 2: Parameterization of the droplet shape by the surface vector R. The deformed surface can be 
parameterized by R'{x, t) = Rq [1 + £(x, t)] N, where N is the unit normal. 

where 

7 = cos -1 (7) 

A. Fluctuation spectrum 

To investigate the influence of hexatic order on spherical droplets, we study deformations about 
the equilibrium configuration. We expand the free energy F d in a small time-dependent displace- 
ment field 8R(x, t), where 

R'(x,t)=R (x) + 5R(x,t) (8) 

is the deformed surface and Ro(x) is given by Eq. ©. For liquid droplets with hexatic order it is 
sufficient to consider purely normal displacement fields, cf. Fig. El Thus, 

8R(x,t) = Ro((x,t)N(x), (9) 

— * — * 

where N(x) = Rq/Rq is the normal vector of the sphere. The dimensionless function ( can be 
expanded in terms of spherical harmonics 

oo I 

(( X ,t)=J2 E r Ut)YUx). (10) 
1=0 m=-l 
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In the absence of defects, the expansion of F d in ( would be straightforward. On the sphere, 
however, one has to deal with 12 discrete disclination charges which produce a small static icosa- 
hedral surface deformation. We initially neglect this discreteness and approximate these defects by 
a smeared distribution of defect "charge". As will be shown in Sect. Ill Bl corrections arising from 
the discrete nature of s(x) are irrelevant for the oscillation frequencies with I < 6. These 
considerations can be made more precise by expanding s(x) in terms of spherical harmonics 

^ oo I 

s(x) = G + ^ S l™Ylm(x), (H) 

1=1 m=-l 

where G is the undeformed Gaussian curvature, G = 1/Rq, and s lm = Yui=i9i^-im{ x i)- To 
obtain the last equation, we have used the representation of the ^-function in terms of spherical 
harmonics, 

oo I 

5(cos 9 - cos 9')6(<p - <pf) = 2 Y *™( 9 '> <f')Yi m (0, <P)- (12) 

1=0 m=-l 

In Eq. © we initially smear out the disclination charge by setting si m ~ for all I > 0. 

With this approximation, the hexatic order has no influence on the droplet shape and the equi- 
librium configuration is a sphere with a radius Rq. However, the presence of hexatic order with 
a nonzero stiffness constant K A nevertheless has an important effect on the fluctuation spectrum 
of a spherical droplet. 

To calculate lo(1), we adopt the treatment of capillary waves on spherical droplets without 



hexatic order IBal and consider the incompressible Navier-Stokes equation for the fluid of the 
droplet 

Pift+ Pi{v-V)v = -Vp + T]V 2 v (13a) 
V-w = 0. (13b) 

Here, p is the hydrostatic pressure in the presence of an interface and pi, r] and v denote the 
density, shear viscosity, and the velocity of the liquid inside the droplet, respectively. 

For droplets the inertial terms in the Navier-Stokes equation (I13al) dominate and effects of 
viscosity are irrelevant. Upon neglecting the non-linear term, the boundary conditions can be 
treated most efficiently by introducing a velocity potential $ with v = V$. Then, integration 
across the droplet interface leads to 

<9$(r) 



<9$(r) 

Pi 



Ot 



Pv 



dt 



= Ap(x), (14) 

r=(R +CRo) + 
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where p v is the vapor density of the surrounding medium. Eq. (TBI) relates the pressure difference 
between the inside and outside of the droplet to the generalized pressure discontinuity Ap(x) 
caused by the shape displacement. Since divv = one has V 2 $ = 0, a Laplace equation with 
solutions of the form IG6I1 

*(r,x, t ) = ( £^« Y '-M(i)' fa -<*d + (15) 
( T,i, m A> m (t)Y^(x) for r > R>(\ + <). 

The displacement field £ and the velocity potential $ are related by the boundary condition that 
the interface velocity must match the fluid velocity, 



RoC = Row: — -rr- 
ot Or 

Hence, the coefficients in Eq. (fT5l) are given by 



(16) 



r=R +(R 



-,2 r lm{t) 



and 

Af m (t) = Rl^- (18) 
Upon setting 

Ap(x) = VA ftm (t)Y, m (a;) = - V) ^fi^Y im (x), (19) 

rCndr, 

l,m l,m u "™ 

(where rj* m denotes the complex conjugate of r /m and the free energy of the deformed droplet), 
one then finds 

^Plmit) =(j + JZpj) Rlhmit). (20) 

Note the influence of hexatic order appears through the term 5F' d /8ri m in Eq. ( fT9l . 

To determine the coefficients pi m in Eq. (fT9l . one has to calculate the variation of the droplet 
free energy. Details about this calculation can be found in Appendix O For small deviations from 
a spherical shape, the deformed droplet has volume V with 

too I \ 

1=1 m =-i / 

(see Eq. (IC6I) ). Thus, the volume constraint V = V appropriate to droplets can be incorporated 
directly by considering only displacements which fulfill (to leading order in the r/ m 's) 

i 

J- \ -* ^ 

^oo 



oo I 



47T 

i=l m=—l 
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With the constraint of fixed volume, the interfacial contribution to the free energy becomes 

/l x ' 

d 2 x^/g~T5~g = F t + -aR 2 £ £ |r im | 2 (/ - 1)(Z + 2), (23) 
1=1 m =-l 

where Fj = ct47t_Rq is the interfacial free energy of the undeformed droplet. Eq. (l23l follows from 
Eqs. (IC5I ) and d2"2"l) . The hexatic free energy is given by 



oo I 



1=5^ EE 



Fir 



Z=l m=—l 



1(1 + 1) 



(24) 



cf. Eq. Co- 
upon setting F' d = F[ + F' h and 

n m = r{Ut)e-*" (,)t , (25) 
and evaluating Eq. fHJt. we find for the fluctuation spectrum (provided I > and p v <C p{) 

k a (i-m +2)- 



to 



PiRl 



/(/-!)(/ + 2) 



1 + 



1(1 + 1) 



(26) 



Eq. d26T) shows that hexatic order only affects undulation modes in a curved geometry: In the flat 
space limit of large Rq and I ^> 1 with k = I/Rq fixed, one has 



LO 2 ~ 



0" + 



i?0 



(27) 



which exhibits explicitely a hexatic correction to the usual capillary wave spectrum. The hexatic 
contribution, however, drops out as R$ — > oo and we recover the result for capillary waves of a 
flat fluid surface [36]. Thus, it is essential to study deformations of a curved geometry to reveal 
the presence of hexatic order. 

In general, the undulation frequency Eq. d26l) depends on the ratio -^A which for hexatic 
order on surfactant-coated water drops or at the surface of supercooled liquid metal droplets is 
Ka/ctRI ~ (£x / 'RoY \E C / 'ctOq) . This ratio grows to become of order E c /aal ~ 0(1) close to a 
continuous hexatic-to-crystal transition, i.e., when £ T « R (see also Sect. I VII) . 



B. Effect of defects on the spectrum 

We now take the deformations associated with a discrete array of 12 disclination defects into 
account, i.e. we consider nonzero si m with / > in Eq. (ITTT) . We first neglect the possibility 
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of disclination motion. Thus, we assume that on the time scale of the characteristic frequency 
(l26l . the disclinations remain in fixed positions which minimize the hexatic free energy of the 
undeformed sphere. Thus, for N d = 12 disclinations at the vertices of an icosahedron, the positions 
Xi in Eq. © are given by Eq. ©. As the sphere is deformed, the hexatic free energy then changes 
as (up to first order in r/ m ) 

S^F h = -K A jrj2 il ~iafi) 2) s ^ (28) 

1=1 m=-l ^ ' 

This follows from Eq. (IC23I) since 5si m is of order r im and the terms of order si m 5si m vanish since 
the defect arrangement on the sphere minimizes F h . To calculate the variation of the interfacial 
contribution the volume constraint has to be included in the free energy. Thus, by considering the 
modified free energy 

F d = F d + [ dVp(x), (29) 



(where p = p ex — pi n is the pressure difference between outside and inside of the droplet), the 
complete first variation becomes 



6Mp d = o J d 2 x^/gj^)5 {l) g(x) + J d 2 x^g\x)p5 w V(x) 

-K A I d 2 x^gjx)^f^ E ^ ~ia + TT 2) a ^im(x) E rl m ,Yl m ,(x), (30) 

J 1=1 m =-l V ' V, m < 

where S^g is given by Eq. (IC4I) and S^'V by Eq. (IC6I) . The shape equation for quasi- spherical 
droplets with hexatic order becomes 

^ oo I 

p. 



+ 2aH(x) = KA E E (/ 1( ]^^ 2) Y lm (x)8 lm . (31) 

Thus, nonzero coefficients si m affect the mean curvature H(x) of the stationary droplet. This 
equation simplifies upon making the ansatz 

oo I 

H(x) = H + 5H(x) = H + E E h ir» Y iM, (32) 

1=1 m=—l 

with 

2H a + p = 0. (33) 
One then finds the extremal equation for the mean curvature, namely, 

2H(*) = ^^t t (34) 
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FIG. 3: Shape of a liquid droplet with hexatic order and 12 disclinations lying on the vertices of an 
icosahedron (surface deformations associated with the defects are exaggerated). 



On the other hand, a deformed sphere has mean curvature (see Eq. (IC15I) ) 

2H'(x) = -J- + yVZ - 1)(Z + 2)r lm Y lm (x). (35) 

l,in 

Comparison of the last two equations then leads immediately to the static surface deformation 
coefficients 

o _ Ka 
rim ~ Slm aR 2 l(l + iy K J 

for / > and rf m = for I = m = in the ground state. Thus, for Ka ^ the defects deform 
the droplet as indicated in Fig. |3] However, nonzero si m have no influence on the frequencies 
u>(l) for < I < 6. Indeed, icosahedral symmetry insures that si m = 0, and hence rf m = 0, 



unless I = 6, 10, 12, ... 113 811 . so corrections of order si m have no influence on the frequencies ui{l) 
for small I. Thus, provided the positions of the disclinations remain fixed on the time scale of 
an undulation (i.e. 5si m = holds in Eq. (IC23I) ) the dispersion relation d2"6l ) remains valid for 
< I < 6. Because disclination motion is catalyzed by absorption and emission of dislocations 
with mean spacing £t, the disclination diffusion constant is D 5 fa (a / '£t) 2 A)> where the D is 
the particle diffusion constant. For a surfactant coated water droplet of radius R = 1mm one 
has (for a = 10" 3 N/m) u(l = 2) ~ 100Hz and D ~ 10" 8 cm 2 /sec. The estimate (a /£r) 2 - 
10 -2 suggests only minor disclination motion during an undulation period and the spectrum is 
unaffected for < I < 6. 
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III. DYNAMICS AND INSTABILITIES IN CYLINDRICAL GEOMETRIES 

Next, we will discuss the influence of hexatic order on the Plateau-Rayleigh instability 
of cylindrical liquid jets coated with a hexatic monolayer. Because the Gaussian curvature of a 
cylinder is zero, defect-free hexatic order is perfectly compatible with this geometry in the absence 
of deformations. However, the hexatic stiffness constant will resist deformations leading to a 
nonzero Gaussian curvature. The complication of defects in the ground state is absent for hexatic 
order on cylinders. 

The surface vector of a cylinder of undeformed radius R is given by 

R(x) = (Rq cos tp, Rq sin (p, z) , (37) 

with x = (ip, z). See again Appendix lAl for the fundamental geometrical quantities of the cylinder 
in terms of R. 

The free energy of liquid cylinders (jets) is given by Eq. (UJ) where now dA = R dipdz. The 
displacement field can again be chosen to be purely normal, i.e., 

R! = R + R (N, (38) 

where now N(ip, z) = (cos ip, sin ip, 0). The displacement field can be expanded in plane waves 

a V ,z,t) = y £r km (t)e lkz e m *, (39) 

k,m 

where k = 2irn/L, with n = 0, ±1, ±2, ... for a cylinder with length L (we assume periodic 
boundary conditions along the axis of the cylinder for simplicity) and 

oo oo 

5> E £ ■ < 40 > 

k,m n=— oo m=— oo 

We neglect the density outside the jet (p v pi), and make an ansatz for the velocity potential 
inside in terms of cylindrical coordinates (r, ip, z), 

*(r,z,<p,t) = J2 A k m (t)e lkz e^I m (kr), (41) 

k,m 

where I m (kr) is the Bessel function of the first kind of imaginary argument 114 lh . Eqs. d 1 4b and 
(fToT) remain valid and yield 
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where I' m (x) = dl m (x)/dx. 

If the cylinder is deformed, its surface area and volume change. Eq. (IC4I) of Appendix ICl 
implies 



A' - A = 27ri?o^ r 00 
and with Eq. (IC6I) one has 

V - V = 2ttR 



k,m 



(43) 



\ fe,m / 



(44) 



where the k = 0, m = term is excluded from the sum Y^ fc . 
Upon choosing 



(45) 



k,m 



the displacement field keeps the volume fixed (corresponding to an incompressible liquid jet) and 
the difference in interfacial free energy between deformed and undeformed cylinder becomes 



F[-F i = nR La £ ' (R 2 k 2 + m 2 - l) \r km \ 



(46) 



k,m 



Note that this energy difference vanishes for k = 0, m = ±1 deformations, corresponding to a 
uniform sideways translation. 

We next calculate the hexatic contribution to the free energy of the deformed cylinder. Eq. 
(ICl 31 ) with background Gaussian curvature G = leads to the Gaussian curvature of the deformed 
state, namely 

G'(ip,z) = J2k 2 rk m e lkz e m *. (47) 



k,m 



Eq. ©, together with the representation of the Green's function of the Laplacian on the cylinder, 



r(x, x') 



Rr 



gifcz ^irrup ^,—ikz' ^—irrvp' 



2ttL k 2 R 2 + m 2 

k,m u 



(48) 



then yields for the hexatic free energy of the deformed cylinder 



Eq. (fl4l now leads to 



W> - 711 K X^W I 2 k4R ° 
*h- Wo «A^\rk m \ k 2 R 2 + m 2- 

k,m u 



n R f mi^M. - „ - 5F d{rkm} 



(49) 



kl^kRo) 



2nR 2 L8rl m 



(50) 
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where 



Upon setting r^r, 



Ap(z,ip) = ^2p km e ikz e 



inup 



(51) 



k,rn 



„0 p —iuj(k,m)t 
km 



one finds the fluctuation spectrum of a liquid cylinder with 



surface hexatic order, namely 

a I m (kR ) 



oj (k } m) 



piRl I m (kR 



Rnk 



R 2 k 2 + m z - 1 + 



k A Rl 



(52) 



oR 2 k 2 R 2 + m 2 

When Ka = 0, Eq. (l52l shows cu 2 < for m = and -Rq/c 2 < 1- The cylinder thus becomes 
unstable if L > 2-kR , leading to the well-known Plateau-Ray leigh instability J3,Q]. However, 
for Ka 7^ the stability of liquid cylinders is enhanced by hexatic order. The 'fastest growing' 
mode kf, which maximizes (— cu 2 ), is now given by 



d_ 

dx 



r 



x 



x=Rok 



X [ 1 — X — 



K 



aR 2 



A 2 
X 



0. 



(53) 



where x = kR . Thus, close to the hexatic-to-liquid transition where Ka/ctRI — 1, one has 
R^k 2 ~ 0.25 compared with R\k 2 ~ 0.48 for Ka = 0. Thus, the characteristic wavelength 
Xf = 2ir/kf of the undulations of the unstable cylinder is significantly stretched by the presence 
of hexatic order. As we shall see in the next section, hexatic order has a similar effect in stabilizing 
multielectron bubbles against fission. 



IV. HEXATIC DYNAMICS AND FISSION IN MULTIELECTRON BUBBLES 



Next, we will discuss multielectron bubbles in liquid 4 He. These bubbles can undergo both a 
freezing transition and a shape instability. Hexatic order affects both the fluctuation spectrum and 
the instability threshold for fission. The free energy of a multielectron bubble F b = F d + F c is that 
of a droplet (cf. Eq. dU)) with an additional Coulomb contribution, i.e., 



a I dA + -K A \ dA D in 3 D'rij 



1 

Ye 



dA / dA 



, Pp)pQ') 
\x — x'\ 



(54) 



Here^jAx) denotes the charge distribution on the surface and e is the dielectric constant of liquid 
4 He 114211 . In an equilibrium fluid, p = eN/AnRl for a sphere with N electrons. 

The assumptions underlying this theoretical description are: (i) The electrons are restricted to 
the two-dimensional manifold given by the liquid-vapor interface. This is justified since density 
functional calculations (cf. e.g. Il43ll and [ 26]) show that the electrons form a thin layer of thickness 
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5 <C R on the surface with 5 ~ 1 — 4nm. (ii) One can neglect the effects on the shape of the 
charged bubbles arising from: (a) applied electric fields which trap or hold the electrons; (b) 
the movement of the bubbles in the system; and (c) gravity. Assumption (a) is justified since 
typical external electric fields are of the order E ~ 3kV/cm I25L while a charged sphere with 
radius R ~ 10/im and iV ~ 10 7 produces an electric field approximately 300 times larger, 
E sp h = f^r? — lMV/cm. Similar arguments apply to (b) since typical drag forces are much weaker 
than the Coulomb forces. Assumption (c) is justified since R is much smaller than the capillary 



length l c of liquid helium, l c = y2a/ pig ~ 0.6mm. Finally, (iii) we assume that the electrons 
can be treated classically, i.e. corrections arising from quantum mechanics can be ignored. For 
electrons on charged bubbles this is justified since at the melting temperature quantum mechanical 
corrections become only relevant at higher densities, n > 10 12 cm~ 2 H25I1 . 

Upon once again neglecting defects by setting si m = 0, I > 0, the stationary solutions of the 
free energy (1541) are spherical bubbles. To determine the equilibrium radius we use the partition 
function of a noninteracting ideal gas to describe the vapor phase. The free energy of a spherical 
bubble is then given by 

F(V, T) = min \k B TN He (log ^%^| - l) + p ex %R z + AttctR 2 + ^ \ , (55) 



R { ° ne V 47rit: 3 /3 7 '^3" ' ■ 2 Re 

where N He is the number of helium atoms, the thermal wavelength and p ex the pressure outside 
of the bubble. Minimization with respect to the first term only yields, of course, ideal gas behavior, 
i.e., 

PinV = N He k B T, (56) 

where pi n is the pressure inside the bubble. However, by taking all contributions into account the 
equilibrium radius R of a charged multielectron bubble is determined by the Laplace equation 

Pin ~ Pex = 2^- - f Ci )2 ■ ( 57 ) 

R 8%Rq£ 

Thus, one obtains (for p in w p ex ) as typical length scale for R the classical Coulomb radius 

Rl ^ £2!. (58) 

Within the approximations described above, the fluctuation spectrum of a multielectron bubble 
with hexatic order can be calculated. Compared with the discussion of Sect.[n]the only difference 
arises from the Coulomb contribution. As derived in Appendix IC 31 one has for p v <^ pi (thus 
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Rc/R 



FIG. 4: Schematic stability diagram for multielectron bubbles with hexatic order. The stability of a spher- 
ical bubble with radius R and charge eN depends on the temperature T. For hexatics in the temperature 
range T m < T < 2$, bubbles with radii R d < R < R' c are stable whereas for T > Ti they are not. The 
enhanced stability is due internal hexatic order. We expect an even larger region of enhanced stability for 
droplets with crystalline order I45L 

neglecting the density inside the bubble) 

(eN) 2 

A Plm = ±—^(1 - l)r lm . (59) 



By combining the last equation with Eqs. (EOT) . d2"5b and (|2"4"1) one now finds (for I > 0) 



w 2 =-^(i- i)(i+i) 

Pl K 



(/ + 2)-4-^ + 



(60) 



aj R2 /(/ + !) 

For = spherical bubbles become unstable to fission if R < R d , i.e. tu 2 (l — l c ) < for 
i?o<i?c/and / C = 2 0Q. 

For Ka 7^ the stability of charged bubbles is enhanced by the hexatic order of the electrons 
on the sphere. Thus, for T m < T < Tj the unstable mode is still l c = 2, but now co 2 (l c = 2) < 
for Rq < R' c with R' c < Rd, cf. Fig. |4j The icosahedral symmetry of the deformed shape with 
Sim 7^ is too high to have an influence on the fission instability which occurs at I = 2, justifying 
our neglect of disclination defects. 

Because the electrons (which determine Ka) are far apart relative to the helium atoms (which 
determine a) KaIoR 2 c1 will be smaller than for droplets of supercooled liquid metals when R ~ 
R d . For helium-bubbles with iV ~ 10 6 one has R d ~ 10/im and we expect (see Sect. IVTb that 
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-^4- ~ 10 4 . Charged metal droplets also undergo the same fission instability. Thus, it might be 
possible to detect the onset of hexatic order by investigating the stability of charged liquid droplets 
in Paul traps 1 19]. 

In the flat space limit (I — > oo, k — I/Rq fixed) one obtains, with n = eN/AitR^, 

i 2 

e 



TL 

p ex uj 2 {k) = ak 3 - Ank 2 — , (61) 



in agreement with 112511 . As Eq. doTt shows, planar He-interfaces are unstable to deformations with 
k < k c = Ann 2 /as. This is the 'dimple' instability discussed in the Introduction which triggers 
the initial creation of multielectron bubbles. 

V. DYNAMICS IN HEXATIC MEMBRANES WITH A SPHERICAL TOPOLOGY 

As a last application we discuss the influence of hexatic order on the fluctuations of a spherical 
vesicle. The free energy F of a hexatic membrane is given by 

F = F b + F g + F h 

= J dA (2H) 2 + k g J dAG + ^K A J dA D^Wn^ (62) 

where we have neglected spontaneous curvature, and k and k g are the mean and Gaussian rigidity, 
respectively. 

The equilibrium shape of enclosed vesicles (provided fluids inside and outside cannot equili- 
brate on experimental times scales) often corresponds to the minimum of the free energy d62l with 
prescribed surface area A and volume V [46] (see also Ref. J^]). The area and volume constraints 



can be implemented by considering the modified free energy 

F = F + F a + F vol = F + J dAa(x) + J dVp{x). (63) 

Often, a constant pressure difference between out- and inside p = p ex ~Pi n = p > —p < and surface 
tension a are taken as Lagrange multipliers to enforce these constraints. Here, we allow for both 
spatially varying surface tension a(x) and pressure p(x) to enforce local incompressibility (cf. Eq. 
d72l below). In many situations, the average values of a and p will depend on the volume captured 
by a vesicle of a given area at the moment of its formation. 

In the following analysis we only consider surface shapes which are topologically equivalent 
to a sphere. Then, the second term of Eq. (1621 can be neglected, since it is a topological invariant. 
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We again temporarily ignore discrete disclination defects and show later that their inclusion does 
not affect characteristic frequencies with I < 6. 

Because k plays a similar role for vesicles as a plays for droplets, hexatic order should lead here 
to similar effects on the spectrum uj{1). For vesicles the dynamical fluctuations take place at very 
small Reynolds numbers Therefore, the convective and the inertial terms can be neglected 
in the Navier-Stokes equation (I13al) for the surrounding bulk fluid. Thus, here we must solve the 
Stokes equation 

Vp = r]V 2 v, (64) 

where r] and v denote now the shear viscosity and the velocity of the bulk fluid, assumed to be 
the same inside and outside the membrane. The undulating membrane imposes the boundary 
condition 



N ■ v(r) 



f=R (l+O 

where R is given by Eq. © and the membrane velocity by 

d((6,<p,t) 



R ( = Rq 



dt 



(65) 



(66) 



The following analysis simplifies somewhat by imposing boundary conditions for the components 
V(iV • v) and V x v instead of vq and v v . Then, 



N ■ V[N ■ v(r)} 
N ■ [V x v(f)) 



r=i?o(l+C) 



-V-(Ro(N), 



0. 



(67) 
(68) 



r=R (l+C) 

where Eq. (l67t follows from the condition div v = 0. 

We proceed now as follows [46]: For a given (deformed) vesicle shape the bending forces on 
its surface are known. These forces must be balanced by the viscous stresses. The induced flow 
field v can be calculated by using Lamb's solution with these boundary conditions. Details about 
this solution have been given by several authors, cf. e.g. Refs. [|46fl . fl49fl . and 15011 . The main 
formulas are summarized in Appendix ID1 

As shown in Appendix iDl the viscous force n = Il n iV + Il t associated with Lamb's solution 
(see Eqs. dHU) and <|D4|)) is given by Eqs. (lDT2l) . (iDlHl and (ITTRIl . By using Eqs. dEB and (iDlOll . 
one then obtains for the normal and tangential force differences 

(69) 



m x)-U> x) 



^ ^ ( 2 / + l)(2/ 2 + 2/-3) v 



1=1 m=—l 



1(1 + 1) 
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and 

00 1 21 + 1 

Uf(x) - U>(x) = J2H r lm (t)I^ V —^—-VY lm (x), (70) 

1=1 m=-l ' 

where fif and fif refer to forces inside and outside of the vesicle respectively. 

The tangential motion of the fluid along the membrane induces lipid flow within the membrane 
itself. However, the membrane stays locally incompressible. This constraint is enforced by the 
local surface tension in Eq. d63t . which balances the tangential component of the stress vector. 
Then, with 

00 1 

a(x,t) = a + ^2 ^2 a irn{t)Yim{x), (71) 

Z=l m=—l 

one finds (for / > 1) 

21 + 1 

vim{t) = ri m (t)Rori - . (72) 

Since any non- stationary membrane configuration exerts a local force on the surrounding fluid 
the shape displacement leads to a generalized pressure discontinuity Ap(x) which balances the 
normal component of the pressure discontinuity between the inside and outside 

U<(x)-U>(x) = Ap(x), (73) 

where Ap(x) can be expanded in spherical harmonics as in Eq. (IT9b . In order to determine the 
coefficients pi m in Eq. (fl9l) now the variation of the vesicle free energy has to be calculated. 
Details are presented in Appendix Id 

The deformed surface is characterized by the Gaussian curvature (see Eq. dC14l l) 

G '( x ) = 4 + 4 E( Z " ^ + 2)n m Y im (x), (74) 

and the mean curvature given by Eq. d35t . The volume change is given by Eq. d2TT) and the change 
in surface area by Eq. dC7l) . which becomes 

too I \ 

2V^r 00 + l^l^ 1 + l ( l + X )/ 2 ) + °( r ?™)- C75) 

1=1 m=-l / 

The volume constraint can again be implemented by considering displacements which fulfill Eq. 
(|22l . Then, upon defining the excess area (relative to a sphere) by 5 A = A — AuR^, the area 
constraint reads 

1 x ' 

o^o Yl Yl l r ^| 2 ^ - + 2) = M = const. (76) 
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Similarly, the area term of Eq. d63l) for the deformed sphere becomes 

1 oo I 

F' a = F a + 2RlY,°i m r* lm + -a RlY, E M 2 (Z - 1)(Z + 2). (77) 

—Ljn __ i=l m=—l 

Finally, it has been shown in 15 ill and H52I1 that the bending energy of a quasi- spherical vesicle is 



given by 



F« b = F b + \kRI I dACCC, (78) 



2 

where F b = 8tik and £ can be evaluated on an undeformed sphere to leading order in (, 



2 

l 



£ = D , D i D 3 D j + -^A- (79) 

Fn 



Upon using the decomposition (flOb of £), one finds 

F 6 ' = F 6 + -kJ] ^ |r lm | 2 Z(Z + l)(Z-l)(/ + 2). (80) 

i=l m=—l 

We now insert Eqs. d77l) . (|24)> . and (l80l) into Eq. (fT9l) (with F' replacing F^) and set again r; m = 
r z°m(^) e_iW( ' / ' )t - Upon equating the pressure discontinuity in (fT9l) to the normal force difference (l69l) 
we find for the fluctuation spectrum of spherical vesicles with hexatic order (for I > 0) 



"(Q = -ijj^(i-l)(i + 2) 



«;Z(Z + 1) + (Tq-Rq + id 



2 , r ^ (/-!)(/ + 2) 



Z(Z + 1) 



(81) 



with 



r(/) = (2i + i){2P + 2i-iy (82) 

Note that cu(l) vanishes for I = 1, corresponding to translations of the vesicle as a whole. 

The eigenfrequencies u(l) explicitly depend on the / = component of the tension a which 
acts as Lagrange-multiplier for the area. Since the value of gq is generally not known, one has to 
use the area constraint (1761 (and the fluctuation-dissipation theorem) to express a in terms of 5 A. 
Which modes pick up the excess area depends on the ratio 7 = ctR^/k. As shown in Ref. [47], 
for floppy membranes (i.e. small 7) each mode contributes equally to the excess area, whereas for 
stiff membranes (i.e. large 7) only the lowest mode (/ = 2) will develop a large amplitude [53]. 

For the purpose of estimating the effect of hexatic order, we proceed in a different way. The 
area constraint becomes much easier to handle if the volume constraint can be neglected. By 
assuming that the vesicle is permeable to both water and larger molecules, the area constraint can 
be incorporated directly (again to leading order in the r i?n 's) by choosing 



2V47T , , 

1=1 m=—l 
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1 

uj ~ — i — 

4:1] 



Then, A' = A and one can set cr = in Eq. (I8TT) . Equivalently, we focus on floppy vesicles, 
formed under conditions such that aRl <C k. 

In the flat space limit of a = 0, large R and I ^> 1 with k = l/R fixed, one has 

Kk 3 + ^k . (84) 

As for the liquid droplets discussed in Sect. HH the hexatic contribution drops out as Ro — > oo 
and we recover the result for undulation modes of a flat fluid bilayer I54L The frequency shift 
(I8TT) depends on the ratio K^j k. However, for large vesicles we expect that — 4k (a universal 
result for flat hexatic membranes at long wavelengths 115 511 ) leading to a frequency enhancement 
in Eq. (I8TT) (with cr = 0) by a factor ~ 13/9 ~ 1.44 for the I = 2 quadrupole mode. Note 
the relatively pronounced effect of hexatic order. In general, we expect the largest change in the 
characteristic frequencies for floppy vesicles with 7 = oR\j k <C 1. 

Similar to droplets, the presence of an icosahedral array of defects in hexatic membranes leads 
to an equilibrium configuration with a deformed surface. For vesicles, the complete first variation 
of the free energy F (see Eq. (1631) ) reads 

5^F =\^\ d 2 x^/gjx) ((2H(x)) 2 5 {1) g(x) + 8H(x)5 {1) H(x)) 



+a / d 2 x^g[x)5 {1) g{x) + / ^x^/gJxjpd^Vix) 



-K A I d 2 x^gjx)±- 2 jr {l ( f^ + 2) s lm Y lm {x) J2 4 m 'n m > (») , (85) 

J 1=1 m =-l \ + > l>, m ' 

where (again) 5^g is given by Eq. d), 5 {1) H by Eq. (IC16Tl and 8^>V by Eq. <|C6l». Since 
integration by parts shows that 

J d 2 x^fgJx)H{x)V 2 C{x) = J d 2 xyfg{x)Q{x)V 2 H{x), (86) 

the equilibrium shape equation for quasi-spherical vesicles with hexatic order becomes now 



p + 2a H(x) - 4H(x)k(H 2 (x) - G(x)) - 2kV 2 H(x) 
00 / 

> ( / — n ( / 4- ' > , . 

(87) 

/ 1 / -f- 1 ; 

1=1 m=—l 



1 ~ • (i-l)(j + 2). 



Upon making again the ansatz (l32l) (where 2Hqo~o + p = 0) and using G = Go + 25H /R , one 
finds for the mean curvature 

= 2^ + I|t + 1) + g fig ~ 2 * s '" Y ''"l :r )' <88) 
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where the static surface deformation coefficients in the ground state are given by 

r?m = Slm KP{l + lf+ A a Rll{l + l) m 

for / > and rf m = for I = m = 0. 

However, icosahedral symmetry again insures that si m = 0, and hence rf m = 0, unless 
I = 6, 10, 12,..., at least for floppy, approximately spherical vesicles. Hence, corrections of order 
Sim have no influence on the frequencies uj{1) for small I. Just as for hexatics on liquid droplets, 
disclination motion is negligible during an undulation period and the dispersion relation (l8lT) re- 
mains valid for < I < 6. This can be seen by using an estimate similar to that of Sect. Ill Bl Here, 
uj{1 — 2) ~ 40Hz for a lfj, vesicle and the disclination diffusion constant is D 5 « (ao/£T) 2 A; P id> 
where D Upid ~ 10 _8 cm 2 /sec and (cio/£t) 2 — 10 -2 . 

Finally, we show that our assumption that the regular part of the bond-angle field relaxes rapidly 
on the time scale of undulation modes is indeed justified for vesicles with hexatic order. As follows 
from Eq. (IB 81) the relaxational dynamics for 9 re9 is described by 

^ m —— = K A V 2 9 re \ (90) 
at 

where i] m is the shear surface viscosity of the membrane. For simplicity, we have neglected in the 
last equation defects which are immobile on the time scale of ripples anyhow. Upon expanding 
9 reg in terms of spherical harmonics 

oo I 

re 9 = J2J2 €f (91) 

1=0 m=-l 

and setting 

€? = 0le~^\ (92) 

one finds 



^(0 = -*^^ + !)- (93) 



K A 

Thus, by comparing with the undulation spectrum uj given by Eq. (I8TT) one obtains with rj m = rjh 
(where ~ 100 — lOOOnm for bilayers) for floppy vesicles and large I 

u e (l) ~ AK A R • 1 ; 

Thus, ojq(1) for h R , justifying our assumption for large vesicles. Similar arguments 

apply to spherical droplets and bubbles and cylindrical geometries. 
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VI. EXPERIMENTAL CONSEQUENCES 

We conclude by discussing the implications of our results for possible experiments. As dis- 
cussed in the Introduction, one of our main motivations was to search for dynamic signatures of 
hexatic order in curved geometries, order which is difficult to detect by other means. As shown 
in Sect. |Vj, the frequency shift due to hexatic order is particularly large for vesicles. The I = 2 
quadrupole relaxation rate can be enhanced by a factor ~ 1.44 due to bond-orientational order and 
should be experimentally detectable by measuring the decay rate of fluctuations by, e.g., video or 
fluorescence microscopy J^fj]. 

For liquid droplets the influence of hexatic order on the fluctuation spectrum is somewhat 
weaker. For the magnitude of the hexatic stiffness one has the estimate |l| 

^ = C(T) (I) 2 ' (95) 

where a is the lattice constant, £r is a translational correlation length which diverges as T ap- 
proaches T m and C{T) is a prefactor. When T > T m one has C{T) ~ 2E c (a/a ) 2 yj], where E c 
is the energy of a dislocation core with diameter a ~ a . On the other hand, for the liquid metal 
droplets a ~ ksT/al and with C(T) ~ ksT one has Ka/<jRq ~ 1 for £t ~ Rq. Thus, hexatic 
order should also have experimentally relevant consequences for the fluctuation spectrum of liquid 
droplets and cylinders when the translational correlation length grows to become comparable to 
the sphere size or cylinder radius. As we saw in Sect.|inl the Plateau-Rayleigh instability will also 
be modified for Ka ^ 0. As the liquid cylinder decays into a chain of droplets the distance be- 
tween the droplets will depend on the ratio Ka/ctRq, as follows from Eq. d52l . For Ka/uRq ~ 1 
one has \f = 2n/kf ~ 12. 6-Rn as typical distance between the droplets in contrast to A/ ~ 9.0-Ro 
appropriate to K A = 

For multielectron bubbles a is determined by the distance b between the helium atoms. Here, 
we expect that a ~ 3/c£T m /&Q, where T m is the melting temperature and bo — 3 A for liquid 
helium. From the experiments by Grimes and Adams |Q] it is known that for the planar case the 
dimensionless ratio T m = tt^P ^ 140 and T m = 0.73K for a charge density n = 10 9 cm~ 2 . 
Thus, for 4 He one has a ~ 3 • lCT 4 J/m 2 at T w T m . For helium bubbles with ~ 10 6 one 
has a = {At[R 2 c1 /N) 1/2 ~ 300A and the critical radius R d ~ 10/im. By using the T = result 
E c ~ 0.1e 2 /ao ILZD one finds K A /aR 2 d ~ 10~ 4 . It is worth mentioning that K A /aR 2 cl ~ l/\/N 
since Ka ~ NE C and E c ~ R^/N 3 / 2 . Thus, the influence of hexatic order will be somewhat 



25 



stronger for bubbles with smaller N. 

Thus, hexatic order has only a weak influence on the stability of multielectron bubbles. But 
with sophisticated experimental setups it still might be possible to measure its effects. A possible 
scenario to study fission could, e.g., make use of the fact that bubbles become compactified as 
they move under the influence of an external field away from the liquid- vapor interface towards 
the positive electrode immersed in the liquid. Thus, if for T > T { multielectron bubbles with an 
initial radius R are stable to fission up to a distance h from the interface then they will be observed 
at a somewhat larger distance h + Sh for T m < T < T ?; . 

It also might be possible to experimentally confirm our theoretical predictions in the regime 
of stable bubbles by measuring the characteristic frequencies of an oscillating bubble which is 
stabilized by the balance of bouancy and electric field forces at a constant height below the surface. 
Our predictions also apply to multiion bubbles which experimentally can be realized by charging 
helium films with ions instead of electrons 112011 . 

Finally, for charged droplets in Paul traps similar effects should be observable by, e.g., trapping 
a droplet and shooting charges on it until it undergoes fission. For this purpose it might be neces- 
sary to use more conventional liquids and to perform the experiments at much higher temperatures. 
Fission with ordered ions on droplets can still be achieved by using ions with high charge Ze since 
T m oc Z 2 e 2 ^/n. For example, for droplets with R ~ 10/im and iV ~ 10 6 the melting transition 
should occur at room temperature for Z = 7. Provided the surface tension of the liquid is large 
enough the droplets will still be stable at these parameter values. 

In summary, we have shown that two-dimensional hexatic order should lead to experimentally 
observable effects in a variety of systems with spherical and cylindrical geometry. Similar effects 
could also occur in related systems such as, e.g., cylindrical vesicles. Here, we expect that the 
laser-induced pearling instability will be modified by the presence of hexatic order 115 6l 15711 . 

Crystalline order will also alter the fluctuation spectrum of spherical and cylindrical droplets 
and membranes [58, The effects discussed here should be even larger if the hexatic 

phase is bypassed and one freezes directly into a two-dimensional solid with shear modulus \i. 
The resulting frequency shifts can be estimated by replacing Ka by fiR 2 , in the formulas above. 
Details will be given in an upcoming publication [45]. 
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FIG. 5: Orthonormal basis vectors eg and e v of the tangential plane of the sphere and basis vectors eL and 
e z of the tangential plane of the cylinder. One has eg = Ri/Rq, e v = R2/R0 sin 9 for the sphere, see Eq. 
©. For the cylinder, e v = R\/Rq and e z = R2, see Eq. d37t . 

APPENDIX A: PARAMETERIZATION OF SHAPE 

Here, we review the elementary differential geometry needed to determine all relevant geomet- 
rical properties, i.e. the first and second fundamental forms, of a surface in terms of its surface 
vector R. For a more extensive review, see [61]. For simplicity, we specialize here to spherical 
and cylindrical surfaces. Deformations of these shapes are treated in Appendix ICl 

For an undeformed sphere with constant radius -Ro one can choose x = (x 1 , x 2 ) = (9, p) with 
polar coordinates 9 and (p. The surface vector is then given by 

R(9, ip) = R (sin 9 cos (p, sin 9 sin ip, cos 9) . (Al) 

The first fundamental form is defined by 

gij = R\ ■ Rj, (A2) 

with the covariant vectors 

R. = R. = dsj R=_. (A3) 

On the sphere the vectors {R\(x), R,2(x)} form an orthogonal (but not orthonormal) basis in the 
tangential plane at the point x, cf. Fig. |5J The contravariant components are given by R k = g kt Ri, 
where is the inverse of g^, i.e. g^gjk = S l k . We use the summation convention throughout. 
With the parameterization of Eq. (IA1I) we have 

7^! = Rq (cos 9 cos p, cos 9 sin tp, — sin 9) , R 2 = Rq (— sin 9 sin <p, sin 9 cos <p, 0) . (A4) 
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The area element is generally given by dA = v fgdx 1 dx 2 , where 

g = det(ftj-). (A5) 
For a sphere the first fundamental form is given by 

M = ° n I (A6) 



i2gsin 2 



and = _Rq sin 



The mean curvature H and the Gaussian curvature G are determined by the second fundamental 
form, which is defined by 

bij = Rij ■ N — —R{ ■ Nj, (A7) 
where N is the (local) unit normal vector to the surface 

NixW)^*^ 1 '^*^ 1 ^, (A8) 

and the second equality of Eq. (IA7I) follows from R4 ■ N = 0. In terms of the second fundamental 
form, and upon raising an index of b^ via the operation g tk b k j = &*, we have 

2H = -b\ and g = det(K) = - = ^fr' j . (A9) 

J 9 det((/y) 

See Fig. For a sphere one has 

— > 

-> R 

N = — = (sin 9 cos <p, sin 9 sin yj, cos 9) (A10) 
and the second fundamental form becomes 

, — Rn . 

M=\ a ■ (aid 

-i? sin 2 fl 

It is easy to check that 6* = —S^/Rq and hence, if = l/i?o an d G = in this case. 

The covariant derivative of a vector with covariant components ai and contravariant compo- 
nents a J is defined by l^l 

Didj = a jyi - afcly and = «^ + a fc r£.. (A12) 

Here, the rf- are the Christoffel symbols of the second kind. They are related to the Christoffel 
symbols of the first kind T ik j by = g kl Tuj, where 

F ik j = ^ (dig k j + djg ik - d k gij) . (A13) 
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d = l/R C 2 = l/R 

d = l/R C 2 = 




FIG. 6: Principal curvatures for spheres and cylinders. In both cases the principal directions of curvature 
are drawn. For a sphere both principal curvatures are identical, i.e. C\ = C2 = 1/R. For a cylinder one has 
C\ = 1/R, C2 = 0. The principal curvatures are the eigenvalues of (—bf). Thus, 2H = —(C\ + C2) and 
G = C1C2. 



Note that the T ikj are symmetric in the first and last index, i.e., 

r ikj = r jki . (Ai4) 

For a sphere with polar coordinates we have 

Yl = C = I*, = = 0, I* = Y e w = - sin cos 9. (A15) 

As can be shown by a direct calculation, the co variant derivatives of g^, g l \ and 5* vanish. Co vari- 
ant derivatives do not commute in general. For scalars one has DiDjf = DjDif if didjf = djdif. 
However, for vectors the commutator of covariant derivatives is given by 

[A,^H = a^- (A16) 
Here, R l kji is the (mixed) Riemann curvature tensor 

R l kji = dAi - diV l kJ + r-r^. - r™r l mi . (Aiv) 

In two spatial dimensions one has the particularly simple relation 

R l kji = 9 km l lm ljiG, (A18) 
where Y j is the antisymmetric contravariant tensor, 

Y j = (5{5i-5i5{)/^g, (A19) 
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FIG. 7: Geometrical interpretation of the Riemann curvature tensor on the sphere. Parallel transport of the 
vector V along the path shown in the figure leads to a rotation (V' — V) 1 = R l k j i V k 5x : >5x l , where 
is the Riemann curvature tensor, defined in Eq. (1A17I >. As a consequence, on a curved surface covariant 
derivatives do not commute. Indeed, [D^DArf = Gm, where G is the Gaussian curvature which in two 
dimensions determines all components of fljL'i- 

with covariant counterpart 7^ = Y^9ik9ji- Note, Eq. (IA18I) shows that the Gaussian curvature 
is determined by the first fundamental form only. This is the content of the famous theorema 
egregium of Gauss. A geometrical interpretation of the Riemann curvature tensor is given in Fig. 
□ 

For a sphere the nonzero components of the Riemann curvature tensor are 

-R212 = "-^221 — — = sin2 0j -R112 = -R121 = — = 1- (A20) 
5n £22 

Covariant derivatives have another useful property: integrals over covariant derivatives of vec- 
tor fields can be integrated by parts, similar to vector fields in flat space. This is the content of the 
divergence theorem: for an arbitrary closed surface M (i.e. a surface with no boundary, dM = 0) 
one has 

/ dADiV 1 = (A21) 

J M 

for any vector v\ Here D(v % = -\di{^Jg)v l ) = divv is the divergence of v on a curved surface. 
Finally, the Laplacian is defined by A = D l Di. It becomes on the sphere 

1 ( 1 d 2 1 d ( d\\ 

Asph = M t^W + ^eoe [ sind d6 ) I ■ (A22) 
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Note that the eigenfunctions of A sph are spherical harmonics 

A-sph^lm = ^2 Y; m . (A23) 

For undeformed cylinders with radius R one can choose x = (ip, z). The surface vector is then 
given by 

R(x) = (Rq cos <p, R sin ip, z) . (A24) 
Here, the first fundamental form is given by 

, s l R l I 

(9ij) = (A25) 



1 



and dA = Rodipdz. 

The second fundamental form becomes 




(6 y ) = . (A26) 



Thus, G = and H = 1/2R . 

Since for cylinders is independent of ip and z all Christoffel symbols vanish, i.e. = 0. 
Thus, di = Di and the Laplacian becomes for this geometry 

1 d 2 d 2 

^> = d ' d ' = T^ + es- <A27) 

The eigenfunctions are plane waves 

A cyl e ikz e imifi = -\(k 2 Rl + m 2 )e ikz e im ^. (A28) 
R 

APPENDIX B: HEXATIC FREE ENERGY 



In this Appendix we derive a useful representation of the hexatic free energy for vesicles, 
droplets and multielectron bubbles. The discussion is kept rather general since it has to be valid 
for both deformed and undeformed spheres and cylinders. The discussion here generalizes earlier 



approaches ||28 



5511 to (arbitrary) curved geometries with defects. See also Ref. H3 511 



Consider, the free energy F h introduced in Sect. HH i.e., 

F h = X -K A j d 2 x^D t n j D t n j . (Bl) 
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We now introduce an orthonormal basis {e a } of the tangential space. (Note that orthonormality is 
not guaranteed for the tangent basis vectors discussed in Appendix 1X1 for example, the "natural" 
basis vectors Rg and for the sphere are merely orthogonal). Then, for an arbitrary vector field 
V one has the contravariant components 

V = V a e a , e a -ep = 8 a p. (B2) 

— * 

Here, Greek indices have been used to distinguish this basis from the {Ri}. The covariant deriva- 
tive of V is now given in terms of 

D i V a = e a -(d l V) } (B3) 

so that 

DiV a = diV a - COipaV 13 , (B4) 

with V a = V^ep ■ e a and 

uJif3 a = ep ■ die a (B5) 

is the affine connection appropriate for this basis [63]. For orthonormal e a one has Uip a = —uJi a p- 
Thus, one can define a covariant vector field A4 by 

Ui a p = e a pA u (B6) 

where e a p = 5\5 2 ^ — 8^81 is the Levi-Civita symbol. 

We now set, n = n a e a . Since rtrii = n a n a = 1 it is possible to introduce an bond-angle field 
6 as angle between n and the local reference frame {eajwith 

n a= i = cos 6, n a=2 = sin 9. (B7) 

Then, di"a a = e a pn^diQ, and it is easy to show that Eq. (IB 1 1) becomes 

F h = l -K A J d 2 x^{x) (D l Q + A') (A6 + A4) . (B8) 

For six-fold bond-oriented order, we identify 6 with 9 + 2n/6, which determines the minimum 
charge of disclination defects. The vector field At is a "vector potential" associated with the 
Gaussian curvature, i.e., 

DiAj - DjAi = diAj - djAi = -G(x)j tj , (B9) 
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see Ref. 116 1H . The last result is a consequence of Eq. (IA18I) and of the representation of the 
curvature tensor in this basis, namely 

R% 31 = d 3 u% - d^ 3 + wjw" - W ?u#. (BIO) 

Note that Eq. (IB9I) allows an explicit construction for the vector potential for an arbitrary closed 
surface with a given Gaussian curvature, that is 

Aj{x) = -g kl ll3 D k {x) J d 2 x'Jg~(F)T g {x,x')G{x'). (Bll) 

Here, T g (x, x') is the Green's function for the Laplacian A g , 

A g T g (x,x>) = 6{x ~ X '\ (B12) 

Note that the Laplacian A g depends explicitly on the metric g, since 

A g f = DWJ = [y/g&f) . (B13) 

We now divide the field into a regular part reff and a singular part Q sm9 which represents 
the contribution of defects. The regular part Q re9 fulfills DiDjQ reg = DjDiQ reg ; we assume that 
this part relaxes rapidly on the time scale of shape deformations. The singular part Q sm s is related 
to the disclination density s(x) by (cf. Eq. ©) 

-"/;,/;/-)-"" = s (x) = ^J2 ^ x ~ (B14) 

V9 i 

where q-i = ±2tc/Q for disclinations in a hexatic. Thus, 

D J (x)Q smg (x)=g kl llj D k (x) J d 2 x' ^fgJ^)T g {x,x')s{x'), (B15) 



where T g (x, x') is the same Green's function as in Eq. dB lift . Note that in the gauge defined by 



Eq. dB 1 II) . one has 

DjyQ*" 8 + D i A i = 0, (B16) 

which is the Euler-Lagrange equation of the functional dB8l> . 

By using Eqs. (IB 1 111 and (IB 151) the hexatic free energy (IB8I) becomes for an arbitrary manifold 
with metric g and Gaussian curvature G(x) 

F h = ~K A j d 2 x^gjx) j dVV^j(G(x) - s(x)) T g (x, x') (G(x') - s(x')) , (B17) 
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where we have integrated by parts and used Eq. (IB 121) . 

In the planar case, the ground state has no defects, i.e. s(x) = 0. Then, using T g (x,x') = 



APPENDIX C: FREE ENERGY OF DEFORMED SPHERES AND CYLINDERS 

As the spherical shape gets displaced its free energy changes. In this Appendix we calculate 
the corresponding contributions arising from (i) the interfacial free energy and the bending energy 
(Sect. IC it . (ii) the hexatic free energy (Sect. IC 2t . and (iii) the Coulomb energy (Sect. lC3t . We 
begin by discussing the geometrical properties associated with shape deformations. 

1. Geometrical properties and variations of the mean and Gaussian curvature 

The displacement of spherical shapes parameterized by Eq. ® leads to new tangent vectors 





(B18) 



% = Ri - RoCb-Rk + RoQN 



(CI) 



Ri(l + C) + Rq&N. 



(C2) 



The first fundamental form then changes according to 



&9ij = 9ij ~~ 9ij — R[ ■ Rj — Ri ■ Rj — — 2-RoC% — Ro( 2 bij + R-lCiCj- 



(C3) 



Correspondingly, the change in the area element is given by 




(C4) 



(C5) 



The volume change is given by 




(C6) 



while the area of the deformed sphere is 




(C7) 
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The spherical harmonic decomposition (fTOb together with Eqs. (IC5I) and (l22l then leads to Eq. 
423J. 

The change in the second fundamental form has only to be known up to first order in (. The 
normal vector of the deformed sphere is given by 

N' = N-RoCRi + 0(( 2 ). (C8) 

Therefore, the second fundamental form changes according to 

(kj)' = R! l3 ■ N' = b l3 + RoDjd + Chj + O(C). (C9) 

With 

(g jk )' = g Jk + 2R (V k + 0(C 2 ), (CIO) 

one finds 

m' = - + RoD*b + o(c 2 ). (en) 

The resulting Gaussian curvature of the deformed surface is given by 

G' = det(6?' ) = - = + 0(( 2 ). (C12) 

9 g + $g 

Therefore, 

G'{x) - G(x) = -2R HG( + ifr/V&yAC* + 0(( 2 ) (C13) 
= -2G(-V 2 ( + 0(( 2 ), (C14) 

which leads, upon expanding G(x) and £(x) in spherical harmonics, to Eq. (173b . 
Similarly, the mean curvature of the deformed surface is given by 

- 2H' = (&!)' = (^)'(^)' = + fy*')^- + Sbtj) + 0(( 2 ). (C15) 

Therefore, 

2H' = 2H - 2R ((2H 2 - G) - i? ^Ci + C(C 2 ), (CI 6) 

which leads via spherical harmonics to Eq. d35l . 

Although not derived here, Eqs. (IClT) . (IC4l) . dCoT) . (IC13I) and (IC16I) actually hold for general 
geometries with mean curvature H(x) and Gaussian curvature G(x). These formulas will be 
needed in the analysis of Sects. HT^BllTTTl and fVl 
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2. Hexatic free energy 



Next, we calculate the change in hexatic free energy. Since Eq. (IB 171) holds for arbitrary 
manifolds one has for the deformed sphere 

F' h = ~K A J d 2 x^g + 5g J d 2 x'^g + 5g 

x (G'(x) - s'(x)) T g+Sg (x, x') (G'(x') - s'(x')) . (C17) 

Here, T g+ s g (x, x') is the inverse Laplacian defined on the deformed sphere, i.e., 

A g+Sg T g+Sg (x,x') = 5( ^l . (C18) 

Furthermore, G' is given by Eq. (IC14I) . As the surface gets displaced the position of the defects 
might change. This can be taken into account by introducing a disturbed defect distribution 

oo I 

I If, 



S'(x) = s(x) + —^ y] yZ ^ S lm^lm{x) 
u l=Q m=-l 

= °° + if ^ s ooY o(x) + Y^[si m + 5s lm ]Yi m (x)J . (C19) 
The coefficient 5s o is determined by disclination "charge conservation", i.e., 

J d 2 xy/g + 5gs'(x) = 4vr. (C20) 

Thus, 

5s 00 = -2r 00 . (C21) 

Since the coefficients 5si m are of order (, s'(x) and G'(x) differ only in order (. Therefore, one 
can replace in Eq. (IC17I) y/g + 5g by y/g and T g+ $ g (x, x') by T g (x, x'). Since for the sphere [35] 



oo I 



Y lm {e>, v ')YUe\ v ") 

1=1 m=—l 

one finally gets 



IV, x") = - E E m l it ' (C22) 

1=1 m=~l ^ ' 



n = \ Ka t f ' r '"'"" 1> "t, 2 i"| S '"' + fe "" )|2 + °(C 2 ^-). <C23) 

1=1 m=~l 

a result needed in Sects. UllandlVl 



2 ^ ^ 1(1 + 1) 

1=1 m=-l 
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3. Coulomb energy 



The force (per unit area) acting on a dielectric (with dielectric constant e) contributed by an 
arbitrary electric field E is given by J27I 

E 2 



P 



(C24) 



Then, in terms of the electrostatic potential if), E = —Vip with V 2 t/> = 0, and one can make the 
ansatz (for r > Ro(l + £)) 

1+1 



if, 



eN 



r 



(C25) 



The liquid-vapor interface has to be a surface of constant potential. Thus, if>(r = R + 
RoC) =const. Upon setting Q m = cr/ m , where c is a constant, one then finds 



_ eN 

Clm Tim 
Kq 



and 



Thus, 



E = —e r 



V 



d_ 

Or 



eN 



if>(r) 

r=R +tR n D 



R? 



I 1 + - l)r /m Y; m J e r . 

\ l,m J 



(C26) 



(C27) 



(C28) 



APPENDIX D: LAMB'S SOLUTION FOR SPHERICAL RIPPLES 



Here, we present the main features of Lamb's solution for obtaining the fluid velocity fields 
from the surface stresses. We follow here closely the presentations in 114611 . 14911 . and H50I1 . where 
more details can be found. 

The inner solution [i.e. for r < _R (1 + C)l °f tne Stokes equation (loH) with viscosity 77 is given 

by 



1=1 ^ 



2 V (l+ l)(2/ + 3) 
with the velocity potential function 

1 



1 + 3 2 / 

r vpi — 



r/(/ + l)(2/ + 3) 



(Dl) 



(pf(r,x, 



m=—l ^ 



r 

Rq 



(D2) 
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and the hydrostatic pressure p < = YliPf{ r i x > w i m 



m=—l 



r 

Ro 



(D3) 



The outer solution can be obtained by performing the replacement I — > — (1 + 1) in the formulas 
above. Thus, 



l=i 



1-2 



l + l 



Here, 



2+1 



and 



pf(r,M) = £ pUt)Yim(x) (^-) 

m=-l ^ ' 

The boundary conditions become 

N ■ v(r) 



i+i 



(D4) 



(D5) 



(D6) 



J I 



I 



4CH0 = ^ i 2^T3) + R^ > YUX) 

l,m 



(D7) 



and 



N ■ V[iV • #(f)] - V • u(r) 



r =J R (l+C) 



(D8) 



By comparing Eqs. (ID7I) and dD8l> with the right hand sides of Eqs. (l65t and (l67b one then obtains 

V< - -r, + 3)0-1) p > _ , (2/ -1)(* + 2) 

Am — '2m'/ ^ > Pirn,— 'Im'l I -\- \ 

and 



I 



21 ' rtm '""" U 2(Z + 1) 
The stress vector associated with this velocity field is given by 

n = rwv + i% 

-* fdv v\ n^.^ . 

= -Np + 7] + -V{r-v). 

\or r } r 



(D10) 



Here, the inner normal component is given by 

n< = ff<(r = (r + CRoD = £ { 2 |^ - 



< 

Im 



I 3 - Al - 3 
(/ + l)(2/ + 3) 



Pim 1* Y, m (a;). 



(Dll) 



(D12) 
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The outer normal component reads 

n> = n<(f = (R +(R r) = £ { 2 ^( / + 2 )(* + l )tfm ~ Pfmj YUX). (D13) 

Finally, the difference between the inside and outside tangential component is given by 

nf - nf = E { 2 jr ((* - 1 + (* + 2VfJ 

Z,m 

+ (7TWT3)»» " W^Tf-) VY "» W ' <D14) 
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